CUDA-Accelerated Geodesic Ray-Tracing for Fiber Tracking

Diffusion Tensor Imaging (DTI) allows to noninvasively measure the diffusion of water in fibrous tissue. By reconstructing the fibers from DTI data using a fiber-tracking algorithm, we can deduce the structure of the tissue. In this paper, we outline an approach to accelerating such a fiber-tracking algorithm using a Graphics Processing Unit (GPU). This algorithm, which is based on the calculation of geodesics, has shown promising results for both synthetic and real data, but is limited in its applicability by its high computational requirements. We present a solution which uses the parallelism offered by modern GPUs, in combination with the CUDA platform by NVIDIA, to significantly reduce the execution time of the fiber-tracking algorithm. Compared to a multithreaded CPU implementation of the same algorithm, our GPU mapping achieves a speedup factor of up to 40 times

A simplified algorithm for inverting higher order diffusion tensors

In Riemannian geometry, a distance function is determined by an inner product on the tangent space. In Riemann-Finsler geometry, this distance function can be determined by a norm. This gives more freedom on the form of the so-called indicatrix or the set of unit vectors. This has some interesting applications, e.g., in medical image analysis, especially in diffusion weighted imaging (DWI). An important application of DWI is in the inference of the local architecture of the tissue, typically consisting of thin elongated structures, such as axons or muscle fibers, by measuring the constrained diffusion of water within the tissue. From high angular resolution diffusion imaging (HARDI) data, one can estimate the diffusion orientation distribution function (dODF), which indicates the relative diffusivity in all directions and can be represented by a spherical polynomial. We express this dODF as an equivalent spherical monomial (higher order tensor) to directly generalize the (second order) diffusion tensor approach. To enable efficient computation of Riemann-Finslerian quantities on diffusion weighted (DW)-images, such as the metric/norm tensor, we present a simple and efficient algorithm to invert even order spherical monomials, which extends the familiar inversion of diffusion tensors, i.e., symmetric matrices.</p

Assessing the feasibility of estimating axon diameter using diffusion models and machine learning

International audienceAxon diameter estimation has been a focus of the diffusion MRI community for the past decade. The main argument has been that while diffusion models always overestimate the true axon diameter, their estimation still correlates with changes in true value. Until now, this remains more as a discussion point. The aim of this paper is to clarify this hypothesis using a recently acquired cat spinal cord data set, where the diffusion MRI signal of both a multi-shell and Ax-Caliber acquisition have been registered with the underlying histology values. We find that the axon diameter as estimated by signal models and AxCaliber does not correlate with their true sizes for axon diameters smaller than 3 µm. On the other hand, we also train a random forest machine learning algorithm to map signal-based features to histology values of axon diameter and volume fraction. The results show that, in this dataset, this approach leads to a more reliable estimation of physically relevant axon diameters than using sophisticated diffusion models

Modified geodesic ray-tracing for diffusion tensor imaging

In this paper we develop a modified ray-tracing algorithm for geodesic tractography in the context of brain Diffusion Tensor Imaging (DTI). Our technique is based on computing multi-valued geodesics connecting two given points and tracking the evolution of adjacent geodesics. In order to do so we introduce a new Riemannian metric given by the adjugate sharpened diffusion tensor, combined with a constraint on the tracts outcome based on the geodesic deviation. We present tractography results, and compare our method with the existing ray-tracing approach and deterministic streamlining. Our preliminary results show an improved performance of modified ray-tracing regarding false positive fibers. We also show experiments on subcortical short association U-fibers, whose reconstruction is well-known to be hard in a DTI setting

Employing visual analytics to aid the design of white matter hyperintensity classifiers

Accurate segmentation of brain white matter hyperintensities (WMHs) is important for prognosis and disease monitoring. To this end,classifiers are often trained – usually,using T1 and FLAIR weighted MR images. Incorporating additional features,derived from diffusion weighted MRI,could improve classification. However,the multitude of diffusion-derived features requires selecting the most adequate. For this,automated feature selection is commonly employed,which can often be sub-optimal. In this work,we propose a different approach,introducing a semi-automated pipeline to select interactively features for WMH classification. The advantage of this solution is the integration of the knowledge and skills of experts in the process. In our pipeline,a Visual Analytics (VA) system is employed,to enable user-driven feature selection. The resulting features are T1,FLAIR,Mean Diffusivity (MD),and Radial Diffusivity (RD) – and secondarily,CS and Fractional Anisotropy (FA). The next step in the pipeline is to train a classifier with these features,and compare its results to a similar classifier,used in previous work with automated feature selection. Finally,VA is employed again,to analyze and understand the classifier performance and results

An innovative geodesic based multi-valued fiber-tracking algorithm for diffusion tensor imaging

We propose a new geodesic based algorithm for fiber tracking in diffusion tensor imaging data. Our algorithm computes the multi-valued solutions from the Euler-Lagrange form of the geodesic equations. Compared to other geodesic based approaches, multi-valued solutions at each grid point are computed rather than just computing the viscosity solution. This allows us to compute fibers in a region with sharp orientation, or when the correct physical solution is not the fiber computed from the first arrival time. Compared to the classical stream-line approach, our method is less sensitive to noise, since the complete tensor is used. We also compare our algorithm with the Hamilton-Jacobi equation (HJ) based approach. We show that in the cases where U-shaped bundles appear, our algorithm can capture the underlying fiber structure while other approaches may fail. The results for synthetic and real data are shown for both methods

Riemann-Finsler multi-valued geodesic tractography for HARDI

We introduce a geodesic based tractography method for High Angular Resolution Diffusion Imaging (HARDI). The concepts used are similar to the ones in geodesic based tractography for Diffusion Tensor Imaging (DTI). In DTI, the inverse of the second-order diffusion tensor is used to define the manifold where the geodesics are traced. HARDI models have been developed to resolve complex fiber populations within a voxel, and higher order tensors (HOT) are possible representations for HARDI data. In our framework, we apply Finsler geometry, which extends Riemannian geometry to a directionally dependent metric. A Finsler metric is defined in terms of HARDI higher order tensors. Furthermore, the Euler-Lagrange geodesic equations are derived based on the Finsler geometry. In contrast to other geodesic based tractography algorithms, the multi-valued numerical solution of the geodesic equations can be obtained. This gives the possibility to capture all geodesics arriving at a single voxel instead of only computing the shortest one. Results are analyzed to show the potential and characteristics of our algorithm

Riemann-Finsler multi-valued geodesic tractography for HARDI

We introduce a geodesic based tractography method for High Angular Resolution Diffusion Imaging (HARDI). The concepts used are similar to the ones in geodesic based tractography for Diffusion Tensor Imaging (DTI). In DTI, the inverse of the second-order diffusion tensor is used to define the manifold where the geodesics are traced. HARDI models have been developed to resolve complex fiber populations within a voxel, and higher order tensors (HOT) are possible representations for HARDI data. In our framework, we apply Finsler geometry, which extends Riemannian geometry to a directionally dependent metric. A Finsler metric is defined in terms of HARDI higher order tensors. Furthermore, the Euler-Lagrange geodesic equations are derived based on the Finsler geometry. In contrast to other geodesic based tractography algorithms, the multi-valued numerical solution of the geodesic equations can be obtained. This gives the possibility to capture all geodesics arriving at a single voxel instead of only computing the shortest one. Results are analyzed to show the potential and characteristics of our algorithm